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Abstract 

A number of intriguing properties emerge upon the formation of the epitaxial interface between 
the insulating oxides LaA103 and SrTiOs. These properties, which include a quasi two-dimensional 
conducting electron gas, low temperature superconductivity, and magnetism, are not present in the 
bulk materials, generating a great deal of interest in the fundamental physics of their origins. While 
it is generally accepted that the novel behavior arises as a result of a combination of electronic and 
atomic reconstructions and growth-induced defects, the complex interplay between these effects 
remains unclear. In this report, we review the progress that has been made towards unraveling the 
complete picture of the SrTiOa/LaAlOs interface, focusing primarily on present ab initio theoretical 
work and its relation to the experimental data. In the process, we highlight some key unresolved 
issues and discuss how they might be addressed by future experimental and theoretical studies. 
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I. INTRODUCTION 



With the advance of techniques to control thin film growth on the atomic scale, the study 
of epitaxial oxide heterostructures is a rapidly developing area of materials science jjfl. Due 
to the ability to produce a well-defined, single-terminated surface [l, 5], oxide interfaces that 
are nearly atomically sharp can now be fabricated. In many cases, the properties of these 
interfaces turn out to be much richer than those of their bulk constituents pHl 1| . Not only 
are these new interface phases of fundamental physical interest, they are also promising 
candidates for novel devices and technology. 

One of the most interesting epitaxial oxide heterostructures to date is the (001) inter- 

n 

face between LaA103 (LAO) and SrTiOa (STO) [7(. Although both materials are con- 
ventional band insulators in bulk form, among the intriguing phenomena observed at the 
SrTiQq /LaAlOs interface is the presence of a high-mobility quasi two-dimensional electron 
gas 0,13]; which emerges when the LaAlO^ film thickness reaches a critical value of 1.6 nm 



(4 LaA10 3 unit cells) [13|, 



■ ll^. This thickness dependence enables external control of the con- 
ductivity of the heterostructure. Reversible control of the metal-insulator transition has been 
demonstrated via an applied electric field 
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151 ] . In addition to these phenomena, at low 



temperature (~ 200 mK), the SrTiOs/LaAlOs interface becomes superconducting [16H18]. 
Furthermore, though bulk SrTiOs and LaAlOa are both non mag netic, experiments suggest 
that the interface may exhibit some type of magnetic ordering [19]. In addition to fundamen- 
tal scientific interests, the LaAlOs/SrTiOs interface is also promising for the development 
of novel applications in nanoscale oxide electronics. For example, devices have recently been 



fabricated 



20j that exploit the fact that the interface conductivity can be "written" and 



'erased" locally with the aid of an atomic force microscope (AFM) tip 



Despite extensive efforts in both theory 23N371] and experiment 
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21 



22]. 



471 ]. the origins of 



these interface properties are not yet completely resolved. In part, this is due to the fact 
that more than one mechanism may play a role in determining the interesting behaviors. 
Also, which mechanism dominates the behavior in a given sample appears to be sensitive 
to the conditions under which the sample is grown and/or how it is processed prior to the 



experimental measurements 
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|. As a result, direct comparison between experiments can 
be complicated. Nevertheless, progress is being made towards understanding the basis of 
the novel interface phenomena, in particular, the origin of the conductivity. 
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To date, three mechanisms have been proposed to account for the emergence of con- 
ductivity at the SrTiOs/LaAlOa interface. The first mechanism is an intrinsic electronic 
reconstruction due to the polar discontinuity at the interface Q]- The driving force 
behind this mechanism can be understood in the following simple picture. Perovskite ox- 
ides have a generic ABO3 structure, where A and B are metal cations. Along the (001) 
direction, this structure consists of alternating AO and BO2 planes. In the simple ionic 
limit, SrTiC>3 consists of charge neutral layers, (SrO)° and (TiC^) , while LaAlOs is com- 
posed of alternating charged layers, (LaO) + and (A102) _ . The deposition of stochiometric 
LaAlOs onto the SrTiC>3 (001) substrate is thus equivalent to the formation of a chain of 
capacitors in series. The potential across each capacitor is additive and, mathematically, 
diverges with the number of capacitors (i.e., the thickness of the LaA103 film). To offset 
the diverging potential, an electronic reconstruction, known as the "polar catastrophe", is 
expected to occur, whereby half an electron (hole) is transferred to the n-type TiC^/LaO 
(p-type SrO/A102) interface so that the internal electric field through the LaA103 film is 
completely compensated. As a consequence of this charge transfer, the interface becomes 
doped, leading to the observed conductivity. 

The polar catastrophe mechanism is believed to be applicable to the SrTiOs/LaAlOs 
interface because of the multivalent nature of the Ti cations, which can exist as either Ti 3+ 
or Ti 4+ . Therefore, it is possible that this type of electronic reconstruction has a lower energy 
barrier and/or ground state energy than the typical atomic reconstructions that occur for 



more traditional semiconductor polar /non-polar interfaces 43. l50i . 

However, soon after the polar catastrophe mechanism was proposed, several groups [401 — 



42| observed that the electronic properties of the interface are significantly affected by the 
presence of oxygen vacancies, a significant number of which are generated in the SrTiC>3 
substrate during the LaAlOs deposition. These results suggest that oxygen vacancies can 
be responsible for the high sheet carrier densities and mobilities measured in experiments. 
Oxygen vacancies are also proposed 2l| to account for the observed insulating-to-metallic 
transition, via a mechanism involving the creation and annhilation of oxygen vacancies on the 
LaAlC>3 surface. While the concentration of oxygen vacancies is suppressed by performing 
the LaA103 film deposition under high oxygen partial pressures and/or post-annealing in 



oxygen, some residual conductance is still observed 
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|42| . suggesting that oxygen vacancies 



may not be the sole source of the interface electron gas. 
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A third possible mechanism, based on the observed intermixing of cations across the 



SrTiOs/LaA103 interface, has also been suggested 
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381 ] . The driving force for cation 



mixing is assumed to be the reduction of the dipole energy at the interface. As such mixing 



resu. 
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ts in the formation of Lai-^Sr^TiOs 
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which is metallic when 0.05 < x < 0.95 



521 ]. it could also lead to interface conductivity. 



Though each of the three mechanisms can explain particular experimental results, there 
is as yet no unified picture to describe the entire phase diagram of the SrTiOs/LaAlOs 
interface and the relations among the different phases. It is not clear whether the whole 
phase diagram has yet been identified, or whether some or all of these three mechanisms are 
adequate to account for the experimental observations. 

The primary goals of this report are to review the progress made towards understanding 
the SrTiOs/LaAlOa interface using theoretical first principles calculations and to provide 
a theoretical perspective on the existing experimental data. We start by reviewing the 
experimental data in Section [Til highlighting different results in the spirit of stimulating 
further studies. Section II III reviews the state-of-the-art first-principles calculations on the 
SrTiOs/LaAlOs interface and discusses their successes and failures in describing the exper- 
imental observations. In Section IIV| we pose what we think are the important unanswered 
questions about the SrTiOs/LaAlOs interface and in many cases try to suggest additional 
experiments and/or theoretical computations. Our conclusions are in Section [V] The in- 
vestigation of SrTiOs/LaAlOs interfaces is a rapidly growing field. Due to space limitations 
and the rapid publication of new results, we apologize to those whose work is not included 
here. 



II. EXPERIMENT 

For (001) oriented growth of perovskites, two types of SrTiOs/LaAlOs interfaces are 
possible, depending on the sequence arrangement. The first has the atomic planes TiC^/LaO 
at the interface and is referred to as the n-type interface; the second is formed by adjacent 
SrO/A102 planes and is called the p-type interface. Experimentally, it is the n-type interface 
that shows the interesting and novel properties mentioned in the introduction, while the p- 
type interface is found to be insulating. Below, we briefly review the experimental results 
on the n-type and p-type interfaces. For a more complete experimental review, especially 
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FIG. 1: Schematics of various types of SrTiOa/LaAlOa interfaces. The light blue denotes SrTi03, 
the orange denotes LaA103, and the light green denotes SrRuO"3. a) A bare n-type interface 
(Ti02/LaO) with AIO2 terminated surface, b) A bare p-type interface (SrO/A102) with LaO 
terminated surface, c) An np-type interface with SrTiOs as the capping layer, d) A buried n-type 
interface with SrRuC>3 as the capping layer. 



concerning the technical details, please refer to Refs. 



53 



54] 



A. The n-type interface 



The simplest n-type interface (which we call the "bare n-type interface") has a TiC>2- 
terminated SrTiOs substrate on which a LaA103 film is grown, creating a TiC^/LaO interface 
and an exposed AlCVterminated LaAlC>3 surface [3, Q] (see Fig. [T^,). In addition, there 
are other variants. Huijben et al. 55] grew capping layers of SrTiC>3 on the LaA103 to form 
a system with both n-type and p-type interfaces (np- type) (see Fig. [It), while Wong et 
al. [56( grew SrRuOs electrodes as the capping layers on the LaA103 to study the electric 
transport perpendicular to the interface (see Fig. [Hi). These geometries create buried re- 
type interfaces. We briefly review and compare the experiments on these n-type interface 
systems. 



1 . Transport 



The transport properties of the n-type SrTiOs/LaAlOs interface depend on the oxygen 
partial pressure, po 2 , at which the interface is grown, as well as whether and what type 



of additional annealing steps are performed prior to transport measurements. The existing 
experimental measurements of the sheet resistance R s , the sheet carrier density n s , and 
the Hall mobility of samples grown at low oxygen partial pressure (~ 10~ 6 mbar) 
without post-oxidization appear to have converged to the following values: R s ~ 10 -2 Q, 
n s ~ 10 16 - 10 17 cm- 2 and fi H ~ 10 4 cm 2 V^V 1 (see Tabled]). 

Experiments measuring the transport properties of samples grown at low oxygen partial 
pressure (10~ 5 — 10~ 6 mbar) and then annealed in O2 reveal a range of results. Ohtomo 



et al. 



7| do not find any significant change in the transport properties before and after 



annealing, whereas Kalabukhov et al. 



401 ] grew samples under similar conditions, but observe 



an increase in R s and decreases in n s and fin (see Table [I]) after annealing. Basletic et al. [5 



401 ] for pre- and post-annealed 



find similar results to those observed by Kalabukhov et al. 
samples. 

There appears to be general agreement that, for samples grown at high oxygen partial 
pressures (~ 10~ 4 mbar), the sheet carrier density is in the range 10 13 — 10 14 cm~ 2 . Some 
discrepancies arise in mobility measurements of the same type of samples: at low temper- 
atures (~ 5 K), i? s ~ 10 4 Q and fin ~ 10 cm 2 V _1 s _1 are obtained by Ohtomo et al. [?| 
and Herranz et al. 41[, while Kalabukhov et al. 40| and Thiel et al. [l6| report R s and hh 



see Tabled]). 
12| predicts that 0.5 electrons 



values two orders of magnitude smaller and larger, respective . 

It is noteworthy that the polar catastrophe mechanism [?], 
per two-dimensional interface unit cell (0.5e/<2g TO , where osto is the lattice constant of 
SrTiOa) at the n-type interface would completely cancel the internal electric field through 
the LaAlOs film. That is equivalent to a sheet carrier density of 3.3 x 10 14 cm -2 , imply- 
ing that the polar catastrophe should be the dominant mechanism determining the sheet 
carrier density in samples grown at high oxygen partial pressures (~ 10 -4 mbar), while 
other mechanisms, such as doping via oxygen vacancies, may be responsible for the higher 
sheet carrier densities of interfaces grown at lower oxygen partial pressures (~ 10~ 6 mbar). 
However, agreement between the predicted and measured sheet carrier densities in samples 
grown at high oxygen partial pressure is not sufficient to rule out the possibility that other 



mechanisms, such as oxygen vacancies on the LaAlOs surface [21|, also play a role. These 
other possibilities may shed some light on the unexplained discrepancies in the measured 
mobilities in these samples. 

We note that samples grown at very high oxygen partial pressures (> 10~ 2 mbar) or 
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TABLE I: Values of sheet resistance R s , sheet carrier density n s , and Hall mobility [in measured 



from transport in 


different studies. The data 


in this table 


are all taken at low temperature 


(T~ 5 


K). 












growth po 2 (nibar / 
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fiH (cm 2 V 4 s 1 ) 
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O A 1U 
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10 13 
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[40] 
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[40] 




during cooling 
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during cooling 











samples annealed at high oxygen partial pressure [4lJ| have been observed to be completely 
insulating. If post-annealing in oxygen fills oxygen vacancies and repairs the oxygen off- 
stoichiometry, then one would expect that with increasing oxygen partial pressure, the sheet 



7 



carrier density would asymptote to ~ 10 14 cm" 2 , the value predicted for polar LaAlOs and 
an ideal ra-type interface as per the polar catastrophe mechanism. We will further discuss 
this situation in Section IIV Al 

Another factor affecting transport is the thickness of the LaA103 films. Samples grown 
under relatively high oxygen partial pressure (~ 10" 4 mbar) are not conducting unless the 
LaAlOs film thickness exceeds a critical value, referred to as the "critical separation" . The 
two studies investigating the thickness dependence of sheet carrier density report critical 
separations of 4 and 6 unit cells 
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42|, respectively. 



2. Spatial extent and orbital character of the electron gas 



The oxygen partial pressure during the growth process affects not only the transport 
properties but also the spatial extent of the electron gas. Two different distributions of 
the electron gas have been observed. For samples grown at low oxygen partial pressure 
(~10~ 6 mbar) and not subjected to further annealing, the electron gas extends over tens to 
hundreds of /zm into the SrTiC>3 substrate and thus has 3D character. These electrons most 
likely originate from extrinsic dopants (oxygen vacancies) throughout the SrTiC>3 substrate. 
On the other hand, samples grown and/or annealed at a higher oxygen partial pressure 
(~ 10~ 4 mbar) have conducting electron gases that are quasi-2D, confined within a few (to 
at most tens of) nanometers of the interface. These results on the thickness of the electron 
gas are based on a number of methods, including solution of the Poisson equation for the 



interface system 
at the interface 



observation of magnetoresistance and Shubnikov-de Hass osci 



direct imaging using contact-tip atomic force microscopy |57|, hard 



lations 



x-ray photoelectron spectroscopy 
the interface electron gas 16 ] . 



45], and modeling of the superconducting properties of 



In addition, experiments employing x-ray absorption spectroscopy [46] have been able 
to furnish information on the orbital character of the electron gas at the n-type interface. 
Specifically, the degeneracy of t2 g states (originating from the Ti atom in SrTiOa) is removed, 
and the lowest states for conducting electrons are of Ti d xy character. 
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FIG. 2: A comparison of sheet carrier density at the n-type interface from different published 
experiments as well as those from first principles theory. The red circles are data from Ref. 55l |. 



the blue squares are data from Ref. [45l | and the green triangles are data from Ref. [13j]. The 



data from Ref. 



55l ] and Ref. 



13j are extracted from transport measurements. The data from Ref. 
45j are obtained using hard x-ray photoelectron spectroscopy. Orange pluses are values directly 
obtained from first principles DFT-LDA calculations. The dashed and solid curves are fits to the 
first principles data where the LaA103 is treated as a continuous medium with field dependent 
dielectric constant (see Sec. IIIIB ll ) The dashed curve is a model using the bulk DFT-LDA band 
gaps for SrTiC>3 and LaAlC>3 whereas the solid curve is corrected for the underestimation of the 
bulk DFT-LDA band gaps. 

3. Thickness dependence of sheet carrier density 



While it is generally agreed upon that there is a minimum required thickness for the 
LaA103 film before mobile carriers and thus conductivity are observed at the n-type interface, 
the details of this behavior raise issues yet to be resolved. As with the transport properties, 
in general, the thickness dependence of n s is sensitive to the sample growth conditions. A 
comparison of the sheet carrier density determined from different available experiments is 
shown in Fig. [2l along with theoretical predictions from our first principles calculations. 
(Similar theoretical curves are reported by Son et al. [34| .) 
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Thiel et al. 



13| show that after annealing, the sheet carrier density for samples with fewer 



than 4 unit cells (u.c.) of LaAlC>3 is below detectable levels. Above 4 u.c, they observe 
a step-like jump in n s , which then remains in the range of 1 x 10 13 — 2 x 10 13 cm~ 2 and 
shows no further discernible thickness dependence. The discontinuous jump in sheet density 
agrees with the polar catastrophe model of the metal-insulator transition, but the value of 
the sheet density is an order of magnitude smaller than the theoretical value of 0.5e/a| TO , 
or 3.3 x 10 14 cm" 2 . 

Huijben et al. 55] grew np-type interface samples with po 2 of 10~ 5 mbar without annealing 
and observe that the samples are conducting for all LaA10 3 thicknesses studied: n s starts out 
small but finite for 1 u.c. of LaA103 and increases with LaA103 thickness until it saturates at 
1.4 x 10 14 cm -2 for >6 u.c. of LaA10 3 . As Fig. [2]shows, here the saturation sheet density is 
closer to the theoretically expected value of 0.5e/as TO , but the carrier density is nonzero and 
significant for all LaA10 3 thicknesses, which is at variance with a metal-insulator transition 
scenario. 

The above results are obtained from transport measurements. Sing et al. (45] employ 
hard x-ray photoelectron spectroscopy to study samples grown under the same conditions 



as Thiel et al. 



13| but find that the carrier density is already non-zero for LaA10 3 films 2 



u.c. thick. In contrast to Ref. [13| but in qualitative agreement with Ref. [55| (despite the 
differing growth conditions), they observe that n s increases with the LaA103 thickness. The 
authors propose that the finite sheet carrier density of samples with 2 u.c. of LaA103 is 
immobile and does not contribute to transport. This raises the interesting and important 
question of whether some or all electrons at the n-type interface participate in transport, 
which we further discuss in Section IHIB 31 



4- Cation disorder and intermixing 

Despite the high level of epitaxy available with modern growth techniques, some research 



groups have o 
gawa et al 



Dserved that significant cation mixing occurs at the n-type interface. Naka- 
12], using atomic-resolution electron energy loss spectroscopy (EELS), find that 
the n-type interface is significantly rougher than the p-type interface. They suggest that 
the finite width of the electron gas (as dictated by basic quantum mechanics) increases the 
electric dipole energy, which in turn is reduced by exchanging Sr with La across the inter- 
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face. Willmott et al. [38| analyze the interface structure using coherent Bragg rod analysis 
(COBRA) and also observe a graded intermixing over 3 unit cells. In addition, they note 
that La-Sr mixing is evident for ~6 A further away from the interface than Al-Ti mixing. 
Since the ionic radii of La and Sr are more than twice as large as those of Al and Ti, the 
observation of the intermixing is interpreted as a natural tendency of the system, rather 
than an effect due to kinetic disruption during growth. Overall, the basic mechanisms that 
are the driving forces for cation disorder require explanation and exploration, as does the 
question of what procedures can effectively reduce such intermixing to create ideal interfaces. 

5. Magnetic properties 

Experimental studies of the magnetic properties of the n-type interface are an active and 
growing area. To date, some intriguing discoveries have been reported by various research 
teams based on using magnetoresistance measurements as an indirect probe of magnetic 



behavior and ordering. Measurements by Brinkman et al. [19j provide the first experimental 
evidence of magnetic behavior. Since then, additional studies 58, suggest that some 
type of magnetic ordering may occur at low temperatures. However, there are still questions 
about whether the system exhibits hysteretic behavior, the sign of the magnetoresistance, 
and whether the magnetoresistance is isotropic or anisotropic. 

Brinkman et al. [19J studied samples grown under a relatively high po 2 ~ 10~ 3 mbar, 
with sheet resistances > 10 4 Q. They observe isotropic magnetoresistance, which they 
ascribe to spin scattering off localized magnetic moments at the interface, an interpretation 
bolstered by their observation of a minimum in the sheet resistance R S (T) with lowering 
temperature. This minimum is reminiscent of the Kondo effect, which stems from the 
interplay between localized magnetic moments and itinerant charges. Furthermore, they find 
that the magnetoresistance at 300 mK is hysteretic, implying the presence of ferromagnetic 
domains. However, the R S (T) minimum is not present for samples grown at lower po 2 
(< 10~ 5 mbar). In contrast, Reyren et al. [16( grew a sample at a low oxygen partial 
pressure of 6 x 10 -5 mbar, which is then annealed under conditions that (presumably) yield 
well oxidized samples. They find positive magnetoresistance, no minimum in R S (T), and 
no hysteresis in magnetoresistance. The two major differences between the experiments of 
Refs. 19] and [3] are the growth po 2 (10~ 3 mbar versus 10~ 5 mbar) and the thickness of the 
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LaAlOs film (26 u.c versus 8 u.c). The former factor is known to affect the sheet resistance 
and mobility, while the latter modifies the sheet carrier density. 
The SrTiOa/LaAlOs samples investigated by van Zalk et al. (58 



are grown under similar 

conditions to those used by Brinkman et al. [19|- This study reports hysteresis below 
300 mK and magnetoresistance oscillations versus external field intensity B at 50 mK. 
Surprisingly, the oscillations are periodic in yB instead of the standard Shubnikov-de Haas 
1/B periodicity. The authors suggest that this stems from the commensurability condition 
of edge states. They also speculate that the magnetoresistance oscillation is related to 
ferromagnetic ordering, and that the quantum Hall effect might be present. 

Shalom et al. 59( studied samples grown at a lower po 2 of 10 -4 — 10~ 5 mbar. In agreement 



with other observations 



3. 



16, 



40| | . the sheet carrier density is ~10 13 cm -2 and the system 
becomes superconducting at 130 mK. In contrast to Refs. [3] and 58], however, the low 
temperature magnetoresistances are highly anisotropic and have different signs for fields 
parallel or perpendicular to the current. They find no hysteresis down to 130 mK, at which 
point their samples become superconducting. Based on these observations, they suggest an 
anti-ferromagnetic ordering at the interface below a Neel temperature of ~35 K. 



B. The p-type interface 

Experimentally, the p-type interface (AlC^/SrO) is less explored than the ra-type interface. 
Nakagawa et al. [12J suggest that the insulating properties of the p-type interface arise from 
the presence of oxygen vacancies. Unlike the n-type inteface, at which the Ti can exist in 
multiple valence states, the p-type interface has no multi-valent element (e.g., O - is very 
unusual in compounds). Hence, they argue that, in order to compensate the internal polar 
field and potential across the LaAlOs film, an atomic reconstruction such as a high density 
of oxygen vacancies is required. 

In a simple ionic picture, an ideal and unreconstructed p-type interface has 0.5 holes per 
interface unit cell to compensate the polar field in the LaA103 film. An oxygen vacancy can 
potentially donate two electrons, and assuming that both electrons are mobile, a density 
of one oxygen vacancy per four two-dimensional (2D) interface unit cells can precisely ac- 
commodate the holes and leave no carriers at the p-type interface. Using atomic resolution 

□ 

EELS and various fitting procedures, Nakagawa et al. [12J infer the presence of 0.32 ± 0.06 
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oxygen vacancies per 2D unit cell in the SrTiOs near the p-type interface, which is three 
times higher than their inferred value of 0.1 oxygen vacancies per 2D unit cell at the n-type 
interface. 

Although these results support the oxygen vacancy picture, they also raise some questions. 
The first issue is whether both electrons on the oxygen vacancy are in fact mobile and hence 
can be donated to the interface. The experimental finding of 0.32 oxygen vacancies per 2D 
unit cell argues that more than 0.25 oxygen vacancies per 2D unit cell are needed. There is 
also theoretical evidence to support this incomplete donation of electrons [25j, a subject we 
return to in Section IIII C 31 

Regardless of the precise number of oxygen vacancies, in the polar catastrophe scenario 
the holes are created by 0.5e/ag TO electrons being transfered away from the p-type interface 
to the other side of the LaA103 film, whether it be the LaAlOs surface or another n-type in- 
terface with a capping layer. Why do these electrons make no contribution to the transport? 
This question also applies to the n-type and np-type interfaces, which similarly should have 
holes on the LaAlOs surface or at the capping p-type interface, but experimentally show 
only electron-like carriers in transport. We discuss this point further in Section TlV CI below. 

In summary, though the n-type interface has interesting properties that are promising 
for potential applications, understanding the p-type interface is equally important from a 
theoretical point of view. 



III. FIRST PRINCIPLES THEORY 



In this section, we focus on the theoretical work on the SrTiOs/LaAlOs interface system. 
Most of the theoretical literature is devoted to elucidating the origin of the conducting elec- 
tron gas, reproducing the insulating-to-metallic transition, and designing new interfaces with 
more interesting properties. Some successes have been achieved in the comparison between 
theory and experiment. However, discrepancies exist and some important questions are un- 
resolved. Although some interesting results have been obtained from model Hamiltonian 
calculations [60j, we concentrate here on first principle calculations. 
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A. Theoretical methodologies 



With few exceptions, most first principles studies of SrTiOs/LaAlOs interfaces to date 
have focused on defect-free and idealized n-type and p-type interfaces. Aside from the 
desire to understand the basic physics and intrinsic properties of high quality interfaces, the 
primary practical reason has been due to computational costs: including imperfections such 
as oxygen vacancies, cation intermixtures, or other atomic reconstructions requires the use 



of relative! 
Refs. 21 



v 



24 



arge 



25 



simulation cells. Theoretical studies reporting on oxygen vacancies include 



32 1 with one work on modeling intermixing |38[. However, understanding 



the ideal interfaces is both a good starting point for gaining theoretical insight and should 
also be relevant to experiments where it is believed that high quality interfaces have been 
achieved. 

Most of the reported calculations use density functional theory (DFT) 6JJ with the local 
density approximation (LDA) [fn) or the generalized gradient approximation (GGA) 62|. A 
subset have employed the LDA+U method to include approximately some effects of localized 



electronic correlations that are missing from standard LDA or GGA DFT calculations (63]. 
Th ls has been done typically by including Hubbari U and exchange J ejections on the 
Ti-d manifold, although Pentecheva and Pickett [24J apply U on the oxygen p orbitals to 
study the possibility of hole polarons at the p-type interface. The U values e mployed and 



the resulting predictions are not unianmous, ranging from 5 eV to 8 eV [24j, |2 7H29I. |64 |. 
In particular, Zhong and Kelly 29( studied the n-type interface properties as a function of 
U ranging from to above 6 eV and find a variety of geometric, electronic, and magnetic 
ground states versus U. 

The simulation cells used fall into three distinct categories, which provide complementary 
information on these interfaces. The most popular approach has been to employ a symmetric 
superlattice: a number of unit cells of SrTi03 are placed next to a number of unit cells 
of LaA103 along (001) and periodic boundary conditions are applied. In this symmetric 
superlattice approach, both the SrTiOs and LaAlOs films are nonstoichiometric: when the 
SrTi03 has an additional Ti02 atomic layer and the LaAlOs has an extra LaO layer, two 
identical n-type interfaces are formed to make a symmetric system (Fig. [3^); extra SrO and 
AIO2 layers create two p-type interfaces (Fig. [3b). The symmetric superlattice allows for 
simulation of a single interface type in isolation, and leads to relatively small simulation cells. 
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FIG. 3: Schematics of simulation cells for various types of interfaces and superlattices for the 
SrTi03/LaA103 system. The shaded green regions denote vacuum; their thickness is merely 
schematic and not indicative of actual thickness used in the simulation. For superlattice simu- 
lations, there is no vacuum in the unit cell. For interface simulations, the vacuum is inserted to 
separate periodic copies and introduce surfaces, a) n-type superlattices, b) n-type interfaces, c) 
p-type superlattices, d) p-type interfaces, e) np-type superlattices and f) np-type interfaces. 

However, due to the imposed symmetry and the lack of stoichiometry, the LaAlOs film is 
nonpolar. This precludes the study of the polar catastrophe and the evolution of the interface 
properties as a function of LaA103 thickness. Due to the nonstoichiometry, one additional 
electron (for double n-type) or one additional hole (for double p-type) is introduced into 
the symmetric supercell, which means essentially 0.5 carriers per 2D interface unit cell are 
doped at each interface in the calculation, i.e., the sheet carrier density is fixed as dictated 
by chemistry. Physically, 0.5 carrie rs p er 2D cell is the required amount to completely 



compensate the LaAlOs polar field 



121 ]. meaning that the symmetric supercell approach 



provides the properties of the interfaces in the thick LaAlOs film limit. 

Experimental growth procedures generate stoichiometric LaA103 and SrTiC>3 films that 
are thin (on the order of nanometers) and often have exposed oxide surfaces or capping layers. 
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Therefore, the next generation of computations employs simulation cells with stoichiometric 
films and regions of vacuum to break periodic boundary conditions along (001), as illustrated 
in Fig. 12b, EH, andET- The presence of a vacuum region in the simulation cell allows for the 
polar potential developed through the LaAlOs film to drop across the vacuum gap instead of 
through the materials themselves. To simulate the stoichiometric interfaces, a slab of SrTiC>3 
is typically used to model the substrate, a LaAlOs film is placed on this slab, and, depending 
on the system being modeled, the LaA103 can terminate at the vacuum or a further capping 
film is added on top of the LaA103. In the first case, the system has a single buried n-type 
(Fig. [3b) or p-type (Fig. Eli) interface with an exposed surface; in the second case, a pair 
of n-type and p-type interfaces coexist on the two sides of the LaA103 film. We refer to 
these system as the n-type, p-type, and np-type interfaces, respectively. (We note that all 
three types have been fabricated in experiments.) These stoichiometric simulations have 
a polar LaA103 film with an internal field going from the (LaO) + to (A102)~ terminating 
atomic layers. Therefore, such simulations can probe the evolution of the polar catastrophe 
and sheet carrier properties with LaA103 thickness. Furthermore, assuming ideal interfaces 
and surfaces, these calculations represent a faithful reproduction of many of the actual 
experiments. Due to the close proximity of the exposed surfaces to the interfaces, the effect 
of surface defects or adsorbates on the interface properties can also be investigated. 



Finally, there is another type of geometry employed in simulations 
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25 



26 



65|, which 

we refer to as np-type superlattices (Fig. Efe)- These superlattices have stoichiometric SrTiC>3 
and LaAlOs but do not include vacuum. The absence of vacuum makes the electronic 
structures of np-type superlattices different from those of rap-type interfaces (Fig. ET), as we 
discuss in more detail in Section IHIB II with more details. 



B. The n-type interface 

In this subsection, we review the state-of-art theoretical results from first principles simu- 
lations of n-type and np-type interfaces, as well as double n-type and np-type superlattices. 
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1. Critical separation 



As discussed above, an insulating-to-metallic transition has been observed at the 
SrTiOa/LaAlOs interface for LaA103 films thicker than 4 u.c. jisj. The polar nature of 
LaAlOa provides an explanation for this behavior, i.e. the polar catastrophe. This tran- 
sition can be understood by considering the SrTiOa/LaAlOs energy band diagram. The 
energy diagram of the n-type interface is shown in Fig. Hb. Drawn schematically on the left 
side are the valence and conduction band edges of the SrTiC>3 substrate, which are composed 
of O p- and Ti d-states, respectively. To the right are the LaAlOs band edges, which are 
composed of O p- and La d-states, respectively. Due to the polar nature of LaAlOa, there is 
an electric field through the stoichiometric LaAlOs film. Consequently, the bands slope up 
linearly inside the LaAlOs film. The upwards slope in Fig. Hb indicates that the positively 
charged (LaO) + layer at the n-type interface attracts electrons. As the figure illustrates, 
the energy gap of the entire n-type interface system is determined by the difference between 
conduction band edge Ti d-states of the SrTi0 3 substrate and the valence band edge O p- 
states at the LaAlOa surface. Each added LaA103 unit cell increases the surface valence 
band edge by approximately a fixed amount, thereby reducing the total energy gap of the 
entire interface system. The smallest thickness of LaA103 that closes the energy gap is 
thus the critical separation at which the system becomes metallic. Fig. shows that the 
np-type interface behaves similarly, although in this case the total energy gap is determined 
by the conduction band edge Ti cf-states of the SrTiC>3 substrate and the valence band edge 
O p-states in the capping SrTiC>3 layer. 

The insulating-to-metallic transition has been directly confirmed by DFT studies on 



stoichiometric interfaces such as those shown in Fig. [3b and [3f |28l . l32H34j . By computing 
local densities of states in each atomic plane, it is straightforward to locate the band edges 
and to recreate energy band diagrams of the style shown in Fig. and [lb. These diagrams 
show precisely the behavior described above. 

While the theoretically predicted and experimentally observed values for the critical sep- 
aration are similar, they are not in quantitative agreement. From the discussion above, one 
can see that the value of the critical separation depends on the band gap of SrTi03. As is 
well known, band gaps are systematically underestimated in DFT calculations. Nevertheless, 
direct comparison with experiment can be achieved by including a manual correction to the 
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FIG. 4: Schematic of energy diagrams for a) rap-type, b) n-type and c) p-type interfaces, respec- 
tively. The shaded parts are filled valence bands. The conduction bands are empty, and their edges 
are denoted by the solid lines. The dashed arrow refers to the direction of charge transfer. 



band gap, which physically amounts to adding 1 or 2 unit cells of LaA103 to the theoretical 
critical separation. In previous work 32], we employ LDA to study the critical separation 
of the n-type and rap-type interfaces. We determine that the increase of the potential with 
each added LaA103 unit cell is ~ 0.7 eV, which corresponds to an average electric field of 
0.19 V/A in the LaA10 3 film prior to the polar catastrophe. The band gap of SrTi03 in 
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our calculations is 1.9 eV. Taking into account the difference between this and the experi- 
mental SrTiOs band gap of 3.2 eV, we predict that the experimental critical separation is 6 
and 5 LaAlOs u.c. for the n-type and np-type interfaces, respectively. Using GGA with a 



similarly underestimated SrTiOs band gap of 2.0 eV, Pentcheva and Pickett 



33j predict a 



critical separation of greater than 6 unit cells for the n-type interface based on a potential 



shift of 0.4 eV per added LaAlOs unit cell. Cen et al. |2l] find similar values as well: us- 
ing GGA calculations with a SrTi03 band gap of 1.8 eV, they find the theoretical critical 
separation of the n-type interface to be 4 u.c. of LaAlOs and predict that the experimental 
critical separation should be larger than this value. Therefore, for ideal n-type interfaces 
and polar LaAlOs films, LDA and GGA-based DFT calculations have reached a consensus 
that the critical thickness is larger than 4 unit cells and is most likely 6 unit cells for the 
n-type interface. 

Comparison to experiment is complicated by the experimental variances discussed in 
the previous section. The theoretical predictions are larger than the experimental value 
of 4 u.c. of Thiel et al. 13] as determined by transport, or even 2 u.c. from hard x-ray 
spectroscopy results of Sing et al. [451] for n-type interfaces with exposed LaAlOs surfaces. 



55], the comparison is better in 



For the buried np-type interfaces studied by Huijben et al. 
the sense that the sheet carrier density saturates for LaAlOs films 6 u.c. thick. However, 
those experiments find finite carrier densities even for 1 u.c. of LaAlOs, in disagreement 
with the metal-insulator transition scenario. 
At the LDA+U level, Lee and Demkov 



28] use a Hubbard correction on the Ti d manifold 
with [7=8.5 eV to shift the conduction band edge of SrTiOs and increase its band gap 
to 3.2 eV. They report a stronger polar field of 0.24 V/A. They simulate three and five 
u.c. thick LaAlOs films and find that the five unit cell case undergoes an metal-insulator 
transition. Therefore, they report a smaller critical separation of 4 or 5 LaAlOs unit cells. 
This value is closer to the experimental one determined from transport measurements on 



the n-type interface with exposed LaAlOs surfaces of Ref. [13] but differs from the other 
two experiments. 

We note that in all the above simulations, the LaAlOs film itself is modeled using standard 
LDA or GGA methodologies. Therefore, prima facie it is not clear why there is a spread of 
predicted polar fields since the polar field is a property of the LaAlOs film itself. To shed 
some light on this issue of variations in the theoretical predictions, we remark that it is 
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FIG. 5: Schematic showing how an LaA10"3 film can be approximated as a series of parallel-plate 
capacitors and the resulting distribution of electric field and potential. The solid lines refer to the 
actual electric field and potential generated by the capacitors. The dashed lines are the averaged 
electric field and potential over each unit cell. 

possible to make an a priori estimate of the polar field in LaAlOs using only bulk properties 
of the material. Denoting the lattice constant of LaAlOs as (XlaOj the LaA103 atomic planes 
along (001) are alternately (LaO) + and (AIO2) - and separated by a distance of olao/2. A 
simple parallel plate model can be made, as illustrated in Fig. EJ The ±e/ag TO surface 
charge on each plane is spread over a SrTiC>3 unit cell area (asTo) 2 and the resulting electric 
fields are screened by the LaAlOs dielectric constant, clao • Since the field only acts in half 
of each unit cell, we have for the potential energy change AV across a unit cell 

Aire 2 



AV 



OLAO 



2 e L Ao(«STo) 2 
or converting to the average internal polar electric field 

2ne 

polar 1 \n 

eLAOlOSTOj 

Using experimental values of osto = 

3.905 A, 

fl LAo = 3.789 A, we have 



(1) 



(2) 



AV 



22.5 eV 



E. 



polar 



5.94 V/A 



(3) 



CLAO e LAO 

Therefore, each calculation of AV or E po i ar , and thus the resulting critical separation, can 
be directly mapped to a value of £lao implicit in the calculation. 
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FIG. 6: a) Electric field dependence of the bulk LaA103 dielectric constant. The solid line is 
a fit to the first principles calculation, using the Landau-Devonshire formalism Eq. (TjTJ) 6(| with 
eo = 40.95 and Eq = 0.15 V/A. b) Strain dependence of the LaA103 dielectric constant. The 
internal electric field is set to be 0.026 V/A (weak enough that the LaAlC>3 remains insulating). 
The c/a ratio is optimized for each a. The strain is measured with respect to the theoretically 
predicted lattice constant of SrTiO"3 (in experiments the LaA103 film is coherently strained to the 
SrTi0 3 substrate). 

In Table [Til we compile the published values of AV or E po i ar and back out 6lao as per 
Eq. (j3J). As expected, larger AV or E po i ar lead to smaller predicted critical separations, as 
well as smaller implied clao The range of clao spans a factor of more than two, which 
raises the question of why the spread is so large. 

The actual value of the dielectric constant of LaA103 depends on the strain state and 
internal electric field in the material, as well as technical choices such as the pseudopotential 
and basis set (a similar dependence study of 6sto is presented in [67|). We present the field 
and strain dependence of clao f° r bulk LaA10 3 in Fig. [6^ and [6)d. The field dependence is 
calculated in a slab geometry because only thin films can accommodate large internal electric 
fields of around 0.2 V/A without becoming metallic. The dependence can be described 
phenomenologically using the Landau-Devonshire formalism [66] (the solid line in Fig. [6^) , 
which yields, approximately: 



e(E) ~ e 1 




-1/3 



(4) 
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TABLE II: Polar properties of LaA103 films and n-type interfaces from first principles calculations. 
Values of the potential change per added LaA103 unit cell A V and the polar field E po i ar are related 
by Eq. ([3]). In each reference, only one of the two values is reported. The same equation then 
provides the dielectric constant 6lao- 



AV (eV) Ep i ar (V/A) Estimated critical separation elao Reference 

0.4 0.1 >6 56 Pentcheva and Pickett [33] 

0.6 0.16 >4 38 Cen et al. [21] 

0.7 0.19 6 32 Chen et al. [32] 

0.8 0.2 5 30 Son et al. [34] 

0.91 0.24 4 to 5 25 Lee and Demkov [28] 



where eo and E are two fitting parameters. 

In addition, £lao is affected by the strain state of the LaAlOs film, which is in epitaxial 
relation to the SrTiOs substrate. The strain dependence versus in-plane lattice constant 



is calculated using the Berry phase method 68j with a weak internal electric field (0.026 



V/A) so that the system remains insulating and the polarization as a response to the ex- 
ternal electric field is in the linear region. The c-axis is optimized for each choice of a-axis 
parameter. The data in Fig. [6] show that €lao is sensitive to the in-plane strain and the 
internal electric field of the LaAlOs film. Both properties depend on technical choices: e.g., 
different choices of basis sets and/or pseudopotentials change the lattice constant of SrTiOs 
by ± 1-2%, which appears to be a minor difference but in fact has a significant effect on 
clao We believe that a large part of the theoretical spread is caused by these sensitivities. 

To summarize, theoretical work to date has shown that the polar catastrophe mechanism 
does indeed predict a critical separation for n-type and np-type interfaces. Furthermore, 
reasonable values of 5-6 unit cells of LaAlOs are found. The well known underestimation of 
band gaps from LDA or GGA means that the actual critical separations must be larger than 
the raw values produced by such approaches. That the predicted critical separations are 
larger than those found in most, but not all, experiments suggests that other mechanisms 
may play a role before or in conjunction with the polar catastrophe. 

Before concluding the discussion of the critial separation, we mention that an insulating- 
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FIG. 7: a) The experimental stoichiometric np-type superlattice and its corresponding energy 
diagram. The electrostatic potential is stair-like, b) The theoretical stoichiometric rap-type su- 
perlattice with periodic boundary conditions and no vacuum region, and its corresponding energy 
diagram. Due to periodic boundary conditions, the potential has a zigzag shape. The unit cell 
for the periodic boundary conditions is illustrated by the red dot-dash line. In both cases, dashed 
lines with arrows show the direction of electron transfer when the polar catastrophe occurs in each 
system. 

to-metallic transition is also recently found in the simulations of np-type superlattices with 
no vacuum regions [65i] . We would like to point out that there is a subtle difference between 
experimental and simulated np-type superlattices (illustrated in Fig. [7]). For an experimental 
np-type superlattice, the films are of finite size, terminate at surfaces, and are not subject to 
periodic boundary conditions. For such systems, a "global" insulating-to-metallic transition 
can occur, which is missing in simulations with enforced periodic boundary conditions. As 
per Fig. [7^,, before any electronic reconstructions takes place, the energy diagram of an 
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experimental np-type superlattice is stair-like: the potential is constant in the nonpolar 
SrTiOs and increases linearly in the polar LaA103 layers. When the superlattice is thick 
enough, the electrons filling the LaA103 valence bands on the surface tunnel into the empty 
conduction bands of SrTiC>3 substrate. This is the "global" insulating-to-metallic transition 
as the transferred electrons and holes generate an additional internal electric field throughout 
the entire superlattice. This additional field counteracts the polar field of LaAlOs film and 
distorts the electrostatic potential, as illustrated in Fig. [7b. Thus the limit of infinitely 
thick experimental np-type superlattices can be simulated by imposing periodic boundary 
conditions imposed on one superlattice unit cell. Separately, if one increases the thickness of 
the LaA103 film in each periodic unit, there will be a "local" insulating-to- metallic transition 
in which electrons in the LaAlOs valence states tunnel into the SrTiC>3 conduction bands 
in the same unit superlattice. We note that the "global" transition (in experiment) or the 
periodic boundary conditions (in computation) generates a counteracting field that polarizes 
the SrTiC>3 and reduces the polar field in the LaAlOs film. This explains the smaller value of 



0.057 V/A 



cells 
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65 [obtained in this approach. In earlier calculations using relatively small unit 
25l |26| , the theoretical np-type superlattices turn out to be insulating at LaAlOs 



thicknesses that make the np-type interface metallic. By employing much thicker films in the 
superlattice unit cell (12 u.c. of SrTiOs and LaAlOs ) an insulating-to- metallic transition 
is found (gsJ, and as expected from the above reasoning the critical separation of np-type 
superlattices is quite a bit larger than that of np-type interfaces. 



2. Orbital character of interface bands 
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30 



32 



34j and experiments |46|, it has been seen that the Ti- 



In both simulations 

derived t 2g degeneracy is split at the n-type interface. Theoretical studies show that the 
lowest occupied bands at the interface reside primarily in the SrTiOs and have a clear Ti 
d xy character. With increasing energy (or equivalently, sheet carrier density), d xz and d yz 
derived states also become partially occupied. Therefore, for all but the very lowest sheet 
carrier densities, multiple interfacial bands, with possibly different orbital character, are 
occupied. 



Detailed band structures for the n-type interface have been presented in Refs. 30 



with some differences. Popovic et al. 30] use a symmetric double n-type superlattice, which 
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FIG. 8: The band structures of an n-type interface and an n-type superlattice from DFT-LDA 
calculations. The solid horizontal red line is the Fermi level, a) n-type interface with 11 u.c. of 
SrTiC>3 and 5 u.c. of LaA10"3 . The thickness of vacuum separating periodic copies is more than 
30 A. The sheet carrier density is 0.1e/a| TO . b) n-type symmetric superlattice with 12.5 u.c. of 
SrTiC>3 (with extra TiO"2 plane) and 4.5 u.c. of LaAlC>3 (and extra LaO plane). As described in 
the text, the sheet carrier density is forced to be 0.5e/a| T Q. 

enforces the presence of precisely 0.5e/a§ TO at each n-type interface (equivalent to what 
would be found for a very thick LaA103 film). Son et al. 34J use a stoichiometric n-type 
interface calculation and present the band structure for 5 u.c. of LaA103, which gives rise 
to 0.2e/ag TO at the interface. The energy difference between the lowest and next lowest 
interface band edges, both of d xy character, is 0.1 eV when 0.2e/ag TO is present at the 
interface, but this value triples to 0.3 eV when 0.5e/a| TO is transferred. In other words, 
there is an evolution of the energy spacings of interface states as a function of electron 
doping. Both works show that the lowest d xy band is localized in the (001) direction and 
thus can be said to have a strong two-dimensional (2D) character. Higher energy interface 
bands extend further into the SrTiC>3 substrate in the (001) direction and are thus more 3D 
in character. 

A convenient way to summarize some of these facts is to show the Fermi surface of the 
interface system in the k x ,k y interface plane. Fig. [8] shows such a diagram. Following the 
above discussion, the lowest energy d xy band has the largest occupancy and thus the largest 
Fermi surface. With increasing energy, the bands have smaller Fermi surfaces. 
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3. Possibility of Anderson localization 



Discrepancies between theoretical predictions and experimental measurements of the 



sheet carrier density 



131 ] have stimulated discussion as to whether all of the transferred 



electrons contribute to the conductivity 30|, |34J. An appealing argument for carriers to have 



different mobilities, first proposed by Popovic et al. 30J, is based on the fact that, formally, 
all electronic states in 2D are Anderson- localized by disorder in the thermodynamic limit. 
Since contributions of localized states to transport are greatly suppressed, the idea is that 
electrons occupying 2D-like states should have a stronger tendency to localize and therefore 
not contribute to transport. 

Popovic et al. argue that the lowest d xy band, which is strongly bound to the interface, 
should therefore be ignored in predictions of the transport carrier density. Furthermore, they 
contend that the higher energy bands with d xz or d yz character should also localize because 
they have large effective masses along the (100) or (010) directions. The only electronic 
states that should contribute to transport are the higher-lying d xy bands. They determine 
the carrier density from these bands to be 8 x 10 13 cm -2 , which is closer to, but still a few 



times 
et al. 



arger than, the carrier density measured by Thiel et al. 13| . The results of Huijben 



55| . however, show a larger value of 1.4 x 10 14 cm 2 . 



Son et al. 



341 ] take this analysis one step further. They exclude the lowest d xy , d xz , and 



d yz bands as well as any d xy bands whose wave function is concentrated within ~ 2 nm 
of the interface, under the assumption that such electrons are easily scattered by interface 
roughness and therefore have low mobility. Only electrons occupying loosely bound or 
unbound d xy bands are counted as contributing to trans por t. They find a sheet carrier 
density that agrees well with the transport data from Ref. 13| (but not that from Ref. 55() 
in that it has a discontinuous jump at a critical thickness and then quickly saturates to a 
value of ~ 2 x 10 13 cm -2 . 

We suggest that a direct experimental test of such localization would be to deliberately 
modify the amount of disorder by growing the same interface system with variable quality. 
Localization would show up as a variation in the transport carrier density. In the absence 
of such data, however, we attempt to estimate theoretically the necessary conditions for the 
likelihood of the Anderson localization scenario at the non-interacting, one-particle level. 
We follow the review by Lee and Ramakrishnan 69]. 
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In 2D, all eigenstates of an infinitely large disordered system are localized: any eigen- 
function ip(r) will behave as |^(r)| ~ exp(— |r — ro|/£), where ro is the site about which 
the state is localized and £ is the localization length. To understand the consequences of 
this for a finite system, the scaling theory of localization |70| describes the evolution of the 
sample conductance, G(L), as the length scale L evolves from microscopic to macroscopic 
dimensions. Localization suppresses the conductance, so that when L > £, 

G(L)~^exp(-L/fl. (5) 

For SrTi0 3 /LaA10 3 interfaces, the theoretical asymptotic sheet density is 0.5e/ag TO , or 
3.3 x 10 14 cm -2 . In order to match transport experiments that give 2 x 10 13 cm -2 , ~90% of 
the electrons must be localized. Their contribution to the conductance must be at least 100 
times smaller than the remaining ~10%. Thus, a conservative bound is exp(L/£) > 100 or 
L/£ > 4.6. In a typical experiment, the contacts used for transport measurements can be 
about 1 mm apart [13], setting another conservative upper bound of L = 1 mm, meaning 



£ < 0.2 mm. 

The scaling theory also provides estimates for £. For the weak scattering scenario, £ is of 
order 

£~Zexp(|M) (6) 

where kp is the Fermi wave vector of the electrons being considered and I is the mean free 
path. (Weak scattering is defined as kpl ^> 1.) The DFT band structures in Fig. [HJ which 



are similar to the results of Refs. 



30 



34j . provide us with kp for each band. From Fig. [8b, we 
can extract out that the lowest d xy band has kp = 0.287r/asTO; while the next higher energy 
d xy band has kp = 0.2l7r/asTO or 2.3 nm _1 and 1.7 nm _1 , respectively. Higher energy and 
less occupied bands have smaller kp. These values, along with Eq. (jBJ), allow us to convert 
the bound on £ to a bound on I. 

If I is approximately the same for all the interface bands, then £ will be larger for the 
lower energy bands, which have the larger kp and will progressively decrease for higher 
energy bands. This is the inverse of the trend needed for localization to hold. Thus, for 
localization to occur, I must be significantly larger for the higher energy bands. Namely, 
since the bands become more delocalized into the SrTiC>3 substrate with increasing energy, 
the interior of the SrTiC>3 substrate would have to be of high quality while the part with a 
few nm of the interface would have to be sufficiently disordered. 
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Following this logic and using the above numerical values in Eq. (jBJ), we find the equivalent 
conditions that / < 3 nm for the lowest band and / < 4 nm for the higher d xy bands. These 
values are small, considering the following standard estimation of mean free path: 



Taking a typical value of Hall mobility (see Tabled]) //# ~ 10 3 cm 2 V _1 s _1 and kp ~ 2 nm -1 , 
then we obtain I ~ 100 nm. Regardless, assuming that a reasonable 2D cross section a for a 
relevant defect is a few nanometers, the typical spacing of such defects in the interface plane 
is ~ \/la and must also be on the order of a few nm or less. 

We briefly consider what defects could produce such small I. One may dismiss surface 
steps or edge dislocations in the SrTiC>3 substrate. The former are typically hundreds of nm 
apart, and experiments indicate that they do not hinder transport, while the latter have a 
small areal density of ~ 10 -8 nm -2 j47|. The presence of extended defects every few nm can 
be excluded, as they would be visible in transmission electron microscopy (TEM) images of 
the interfaces. This means that point-like defects in the immediate vicinity of the interface 
are the most realistic possibility, and candidates include cation intermixtures and oxygen 
vacancies. 



There is experimental evidence of La-Sr intermixture across the n-type interface 12 



on a scale of 1-2 nm. Son et al. 



34j use this "interface roughness" as the basis for excluding 
states localized within ~2 nm of the interface. However, this roughness is in the direction 
normal to the interface, whereas the relevant spacing is that of defects in the interface 
plane. For intermixtures to be responsible, La-Sr mixing must be rough over the scale of 
a few nm in the interface plane itself. Whether that occurs remains to be determined in 
future experiments and theoretical calculations. 

Regarding oxygen vacancies, Nakagawa et al. infer a small density of oxygen vacan- 
cies (0.10±0.04 vacancies per 2D unit cell) at n-type interfaces. This defect density is of 
the right order of magnitude and might be considered a viable candidate for the disorder 
mechanism. Whether these oxygen vacancies have the necessary cross sections for scattering 
is a separate question. More importantly, oxygen vacancies at this density can in principle 
donate 0.20e/cig TO (1.3 x 10 14 cm~ 2 ) to the interface, which is larger than the observed 
transport density. Even if the donation is imperfect, the number of electrons being added to 
the already electron-populated interface would be significant and these electrons would oc- 
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cupy high energy, loosely bound states. They would greatly increase the number of electrons 
contributing to transport. 

In summary, the Anderson localization scenario provides an attractive, yet still unsub- 
stantiated, picture to explain the low sheet carrier densities observed in the transport ex- 
periments of Thiel et al. 47||. A more quantitative analysis of the situation leads us to the 
following conclusions: i) to induce the required localization, the mean free paths of states 
strongly bound to the interface must be smaller than those that are loosely bound; ii) the 
mean free path in the interface plane for the strongly bound states could be at most a 
few nm; and iii) point-like defects are likely scattering centers, with cation intermixtures 
or oxygen vacancies being possible candidates. While theoretically feasible, the localization 
scenario requires that a number of factors work hand in hand. 

We end by noting that this localization scenario jsoj is proposed to explain why theoretical 
values of sheet densities are higher than those found in the experiments of Thiel et al. 13 ]. 
On other hand, the sheet densities measured by Huijben et al. [55] are larger in magnitude, 
saturate to a value only a factor of two smaller than the theoretical value of 0.5e/ag TO , 
and show a monotonic increase with LaA103 thickness, all in much closer agreement with 
existing theoretical findings. 



4- Spatial extent of the electron gas 

As discussed in the review of experiments above, the experimental consensus for the 
thickness of the electron gas at the n-type interface is that it is confined to a length scale 



on the order o 
mbar) 



16 



45. 



nanometers in the samples grown at high oxygen partial pressure (~ 10~ 4 
For samples grown at low oxygen partial pressure (~ 10 -6 mbar), the 
conduction electron gas extends hundreds of /im into the SrTi0 3 substrate {43], ||37]. In the 
latter case, the electrons may be due to extrinsic sources (i.e., oxygen vacancies) and thus 
not determined by the intrinsic properties of the interface. The main contribution of ab initio 
theory in this matter is to provide detailed information about the shape and distribution of 
the electron gas in the intrinsic limit. 

Janicka et al. [3jJ use double n-type symmetric supercells and find that the charge density 
profile of the electron gas follows an exponential decay, e~ z l & (where the first TiC>2 layer is 
at z = and the SrTiC>3 substrate is on the z > side). Fitting to DFT data, they obtain 
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8 ~ 1 nm. Due to the imposed symmetry, the geometry used by Janicka et al. 3l( contains 



five SrTiOa unit cells of substrate (~ 2 nm). In our work [32j, we use a stoichiometric n-type 
interface in a supercell with 11 u.c. of SrTiOs, and we find that the conduction electron 
density profile decays rapidly for the first few unit cells away from the interface, but has a 
longer, non-exponential tail further into the substrate. Son et al. 34J have performed the 
largest simulations to date with 15 to 30 u.c. of SrTiOs and also find exponential decay 
close to the interface that turns into an algebraic decay further into the substrate (see Fig. 
5 of that work). Specifically, their calculations show that within 4 nm of the interface, the 
charge density decays exponentially with 5 = 1.84 nm. Further into the SrTi03 substrate, 
the charge density profile is found to decay more slowly with the approximate form b/ (z—Zq). 



5. Binding mechanism to the interface 

The confinement of the electron gases at interfaces is usually attributed to band bending, 



ail 



find 



a well-known concept in metal-oxide-semiconductor (MOS) interfaces. Janicka et al. 
that the confinement of the quasi 2D electron gas at the n-type interface of SrTiOs/LaAlOs 
is not solely due to band bending. They suggest that the formation of metal-induced-gap- 
states (MIGS) due to the presence of the conduction electrons at the interface plays an 
important role, allowing the electron gas to extend further into the SrTiOs substrate than 
simple band bending would allow. 

An important question is why the electrons are bound to the interface in the first place. 
One answer is that, after the polar catastrophe takes place, the electrons at the n-type 
interface feel the electrostatic attraction from the holes on the other side of the LaAlOs film, 
and this attraction binds them to the interface. While this must play a role, it is not sufficient 
for the following two reasons. First, for the double n-type symmetric supercells, which have 
no polar fields and where there are no holes present to create the electrostatic attraction, 
DFT simulations still find a bound electron gas at the interface. This has already been 



reported in Refs. [30|, |3lJ. In fact, for symmetric supercells with no polar fields, one would 
expect the electrons to occupy the lower-energy conduction band edge of the two materials 
(SrTiOs for this pair of materials) and to spread throughout the substrate. Second, holes at 
the p-type interface should feel the same attractive force, but instead are found to delocalize 
into the SrTiOs substrate. We report this effect for stoichiometric p-type interfaces [32] and 
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have verified that this holds for double p-type symmetric superlattices. Thus, the attraction 
of the opposite carriers across the LaAlOs film is not sufficient to strongly bind them to 
the interface. Furthermore, there is an asymmetry between n-type and p-type interfaces in 
terms of binding carriers that is intrinsic to the interface structure itself. 

We have proposed an explanation for this asymmetry [32) . In both SrTiOa and LaA103, 
the conduction band edge stems from the transition metal d orbitals. However, in SrTiOa 
the Ti atoms are on the B site while in LaAlOs the La atoms are on the A site (in the 
standard ABO3 notation). Thus in the bulk phase of SrTiOa, the nearest Ti atoms are one 
lattice constant apart, but at the n-type TiC^/LaO interface the Ti-La distance is y/3/2 
times one lattice constant. This proximity, combined with the large spatial extent of the La 
d orbitals, enhances the tight-binding hopping element between the two transition metals 
to the extent that it becomes about 100 times larger than the hopping element between 
Ti d orbitals in bulk SrTiC>3. This large hopping matrix element creates a pair of bonding 
and antibonding states composed of superpositions of the Ti and La d states right at the 
interface: the low-energy bonding state is primarily of Ti d character and the anti-bonding 
one of La d character. Compared to the other Ti atoms in the SrTiC>3 substrate, the Ti at 
the interface are special, as they have a lower energy and are thus the most favorable sites 
for electrons to accumulate when the metal-insulator transition takes place. (As detailed 
further below, the p-type interface has no enhanced hopping, and thus has no such special 
binding mechanism.) 



In Ref. 



32], we verify this explanation by mechanically moving the Ti at the interface 



to elongate (weaken) or shorten (strengthen) the Ti-La bond. As expected, when the Ti-La 
bond is strengthened, electrons more strongly localize towards the interface, and when the 
bond is weakened, electrons move farther into the substrate. As an alternate check, we 
compute the on-site energies of the transition metal d xy orbitals from first principles in the 
interface region for a 4 u.c. thick LaA103 film that has undergone the polar catastrophe 
(within LDA). The results are shown in Fig. [9j The on-site energy for atom i is the diagonal 
matrix element (d xy \H\d xy ), where H is the Kohn-Sham Hamiltonian and d xy is the atomic 
orbital at that site. As expected, the La d xy is at higher energy than the Ti d xy due to the 
conduction band offset. The salient point of the figure is that the on-site energy of the Ti 
atom right at the interface is not the lowest one, but simulations show that the electron 
density is highest on the first interfacial Ti atom at the interface. The resolution to the 
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FIG. 9: Green squares show the on-site tight-binding matrix elements at the n-type interface 
computed from first principles for a simulation with 11 u.c. of SrTiOs and 4 u.c. of LaA103 
(the insulating-to-metallic transition occurs at this thickness in the DFT simulation). The on-site 
energies are for the Ti d xy and La d xy orbitals. The La atom is at site 0, and the Ti atoms are at 
sites 1-6. The n-type interface is located between sites and 1. Note that the first Ti atom (at site 
1) does not have the lowest on-site energy, but has the highest electron density. The red triangles 
show the effect of including the large hopping element between Ti and La at the interface. The La 
level is pushed up and the Ti level is pushed down, so inclusion of the hopping element now makes 
the first Ti the lowest energy site for electron accumulation. 

puzzle is to include the large Ti-La hopping element of 0.7 eV. As indicated in Fig. [9l 
this pushes up the La level and lowers the Ti level at the interface. With this hopping 
included, the Ti at the interface is now the most favorable site for binding electrons and 
should accumulate the largest number of electrons. 

These findings imply that the spatial extent of the electron gas is determined in large 
part by the Ti-La coupling at the interface. In particular, one does not expect to find 
an equally strongly bound electron gas when LaAlC>3 is replaced by another perovskite in 
which the transition metal is either on the B site or has more localized atomic orbitals. 
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An experimental test of this prediction is straightforward in principle and requires epitaxial 
growth of other oxides on SrTi03 and the same types of measurements already performed for 
SrTiOa/LaAlOs interfaces. If this physical picture is correct, one can consider a number of 
possibilities for engineering the spatial profile and properties of the electron gas. Chemical 
substitutions right at the interface will change both the spatial extent of atomic orbitals 
as well as the interatomic spacings and thus modify the tight-binding hopping matrix ele- 
ments. Alternately, mechanical motion of the transition metal atoms might be achieved by 
using piezeoelectric or ferroelectric substrates that will in turn change the hopping elements 
through direct modifications of inter-atomic distances. 



6. Magnetic ordering 



As discussed above, experiments have shown intriguing magnetic properties at the 
SrTiOs/LaAlOs interface. Several groups have attempted to tackle th e q uestion of mag- 



netic ordering at the interface theoretically with the LDA+U approach 
Pentcheva and Pickett studied both unrelaxed 



h 24, 3, 



24J and relaxed j27| n-type interfaces 



29 



64]. 



using a c(2 x 2) interface unit cell with U = 8 eV and J = 1 eV. They report that atomic 
relaxation does not play an important role in determining the magnetic states and their 
energetic orderings. The ground state in both cases is a ferromagnetic insulating state with 
disproportionated charge ordering in a checkerboard pattern in the plane of the interface. 
Perpendicular to the interface, they find that the excess 0.5 electrons per 2D cell are localized 



and reside in a single atomic plane right at the interface. Without the Hubbard U [27|], they 
find that the ferromagnetic state and charge ordering are suppressed, that the interface is 
metallic, and that the excess electrons are delocalized over several SrTiOs unit cells. 

Kelly et al. {^J use the rotationally invariant LDA+U method [ti| with U = 5 and 
J = 0.64 eV. In contrast to the results of Pentcheva and Pickett, they find that the atomic 
geometry, in particular a GdFeOs-type distortion, is crucial in reducing the bandwidth of 
occupied Ti d states and in stabilizing magnetic orderings. Their ground state is a p(2 x 2) 
anti-ferromagnetic insulating state with charge ordering. However, they point out that 
by imposing ferromagnetic ordering, they obtain a c(2 x 2) insulating state with charge 
ordering, similar to the finding of Pentcheva and Pickett, which is 10 meV higher than their 
anti-ferromagnetic ground state. However, this energy difference is small, close to the limit 
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of their calculation accuracy 29|. 



Janicka et al. 64| use a (1 x 1) interface unit cell and obtain different results. They find 
that when the excess electrons are forced to be confined to within 1.5 unit cells of SrTi03 at 
the interface, the ground state is a metallic ferromagnet, even in the absence of a Hubbard U 
correction. By increasing the SrTi03 thickness, they observe a decreasing magnetic moment 
and disappearance of the ferromagnetic state for 4.5 unit cells of SrTiOs. Upon applying 
U =5 eV, they find that ferromagnetism is stablized in the ground state for all the thicknesses 
(up to 4.5 SrTiOs unit cells). However, their interfaces remain metallic, in contrast to the 
above results. One explanation is that a (1 x 1) unit cell is too small to allow for charge 
ordering — all Ti sites are forced to be equivalent, the eigenstates must be partially filled 
extended Bloch states, and one thus expects to find metallic behavior. 

Charge ordering is an important ingredient in determining the correct magnetic ground 
state of the n-type interface, and the studies above demonstrate that accurate simulations 
require supercells with sufficient lateral extent to allow for this phenomenon. We end this 
section by noting that, while the LDA+U approaches to date have found magnetic states at 
the interface, the calculations with larger unit cells find insulating ground states, in contrast 
to experiments that report mobile electrons even at the lowest of temperatures. 



C. The p-type interface 

A number of first principles calculations have attempted to shed light on the insulating 
nature of the p-type interface as well as its general physical properties. For the insulating 
behavior, both intrinsic and extrinsic mechanisms have been investigated. Concomitant 
with the smaller number of experimental studies of this interface, there has also been less 
theoretical work on this system. 



1. Intrinsic p-type interface: DFT results 



Theoretical calculations that have used DFT within the standard LDA or GGA approach 



report a metallic state at the p-type interface 



24 



25 



321 ] . The holes at the interface reside 



on the perovskite oxygen sublattice, their electronic states have essentially pure oxygen p 
character as expected for the valence band edge of a perovskite, and they occupy the top 
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portion of the valence band, creating a Fermi surface. At this level of description, the 
metallic character is clearly at odds with experiments. 

As mentioned above, the holes at the p-type interface are not strongly bound to the 
interface but instead delocalize into the SrTiOs substrate. We have found this to hold for 



both stoichiometric p-type interfaces [32] as well as large double p-type symmetric supercells. 
Unlike the n-type interface, where binding is enhanced by the unusual Ti-La tight-binding 
hopping element, there appears to be no such binding force present for the p-type interface. 
The oxygen sublattice is continuous and unmodified across the p-type interface. Therefore, 
one does not expect any special hopping matrix element to be present, and calculations 
report that the oxygen-oxygen hopping elements are essentially constant and unchanged 
across the interface. The polar field in the LaA103 will drive holes out of the LaAlOs film 
and into the SrTiOs film at the p-type interface, but once in the SrTiOs, the holes do not 
have any reason to be bound to the interface. 



In our previous work [32], we investigate supercells with up to 11 u.c. of SrTiOs, and 
find that the holes are delocalized within the SrTiOs substrate for such thicknesses. The 
holes do feel the attractive electrostatic force from the electrons on the other side of the 
LaAlOs film. Therefore, one expects that the holes should be bound to the interface for 
thicker SrTiOs films. What the calculations show is that the length scale for the hole 
binding is larger than 11 SrTiOs unit cells, meaning that the binding energies will be weaker 
as well (when compared to the electrons at the n-type interface). At the DFT level, one 
expects an unbound or very weakly bound hole gas for the p-type interface. However, the 
DFT prediction of a metallic hole gas will still hold for such an ideal interface. For a more 
complete picture of how insulating behavior might arise in this picture, see the section [Til C 31 
on oxygen vacancies further below. 

2. Intrinsic effects at the p-type interface: the hole polaron 
Pentcheva and Pickett 



24( investigate the possibility that strong correlation effects can 
induce a self-trapped hole polaron at the p-type interface. To this end, they use an LDA+U 
approach with U—7 eV on the oxygen p-states. They study a symmetric double p-type 
superlattice with the atoms fixed in the ideal perovskite positions and a p(2 x 2) interface unit 
cell. They consider both ferromagnetic and anti-ferromagnetic orderings. The ferromagnetic 



35 



state they find exhibits disproportionated charge ordering and is half-metallic, at variance 
with the insulating behavior observed in experiments. They find the antiferromagnetic 
state to be charge ordered and insulating, with a very small energy gap of 50 meV. The 
half- metallic ferromagnetic state is energetically favored by 0.15 eV over the insulating anti- 
ferromagnetic state. 

These results provide some evidence that correlation effects missing in LDA or GGA 
approaches can create interface electronic structures that differ from the simple metallic 
prediction of the LDA/GGA. With increasing U, charge-ordered localized states become 
energetically favored due to the increasing cost of double occupancy on an atomic site. 
Therefore, one expects a transition from the delocalized, metallic bands of the LDA/GGA 
prediction to localized states upon increasing U. Open questions are how sensitive these 
results are to the particular U value chosen and what the "correct" value of U may be. 
Separately, the above findings are based on unrelaxed atomic structures, and it is unclear 
whether atomic relaxation would favor antiferromagnetism. 



3. Oxygen vacancies at the p-type interface 

Extrinsic effects such as those from defects or impurities may also help explain the in- 
sulating behavior of the p-type interface. Oxygen vacancies are one of the most prevalent 
defects in the SrTiOs substrate, and thus a number of research groups have focused on 
studying their interactions with the p-type interface. That oxygen vacancies can cause the 



12| 



insulating behavior was first suggested by Nakagawa et al. 

Pentcheva and Pickett 24] consider unrelaxed p-type superlattices with 25% of oxygen 
vacancies (i.e., one oxygen vacancy per four 2D unit cells) in either the SrO or AIO2 atomic 
planes at the SrO/A102 p-type interface and find that the ground state is insulating and 
nonmagnetic in both cases. However, the Fermi level is found to be located in a narrow dip 
of the density of states, and the energy gap is too small to distinguish. (When examining 
oxygen vacancies, Pentcheva and Pickett do not apply the Hubbard U that they use when 
investigating the hole polaron.) 



Park et al. [25] studied both unrelaxed and relaxed double p-type symmetric superlattices 
without Hubbard U. Upon allowing relaxation and including oxygen vacancies in the SrO 
plane at the interface, which favors oxygen vacancy formation over the AIO2 layer, they find 
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that the p-type superlattice with 25% of oxygen vacancies remains metallic, in contrast to the 
prediction of the simple ionic limit. However, increasing the oxygen vacancy concentration 
to 50% results in insulating behavior. They suggest that strong hybridization between Ti 
d-orbitals and O p-orbitals close to the Fermi level makes it difficult for the two electrons of 
an oxygen vacancy to easily compensate the holes. This result is in qualitative agreement 
with the experimental finding that 32% oxygen vacancies per 2D unit cell (and not 25%) 



are found at the insulating p-type interface 



12|. 



The above two studies employ symmetric superlattices and consider oxygen vacancies 
directly at the p-type interface. Our work on stoichiometric p-type supercells with a p(2 x 2) 



interface unit cell 



32| investigates the energetics of oxygen vacancy formation in the atomic 



planes at and away from the p-type interface. We find that, for a 25% oxygen vacancy 
concentration, the formation energy decreases as the oxygen vacancies move away from the 
p-type interface. The convergence to the bulk SrTiC>3 oxygen vacancy formation energy is 
achieved at a distance of three unit cells from the interface. These results show that the 
p-type interface repels oxygen vacancies, which instead prefer to stay in the bulk of the 
SrTiC>3 substrate. 

Since both holes and vacancies are repelled by the interface, these results suggest that the 
interface region is not only defect free but also free of holes. Therefore, the interface itself 
will be insulating. The remaining issue is to understand why the holes do not contribute to 
transport. Based on the picture above, we propose the following qualitative possibility: the 
SrTiC>3 substrate will inevitably have some oxygen vacancies, either due to growth conditions 
or thermal fluctuations. As long as the SrTiC>3 substrate is thick enough to contain at least 
one oxygen vacancy per four 2D unit cell (a very low 3D concentration for a thick substrate), 
the holes can become trapped (i.e. annihilated) by the two electrons of the oxygen vacancy 
and thus will not contribute to transport. This hypothesis could be tested by large scale 
simulations that contain p-type interfaces, oxygen vacancies, and relatively thick SrTiC>3 
slabs. 



IV. OUTSTANDING PUZZLES 

Despite experimental and theoretical efforts, the physical origins of some aspects of the 
SrTiOs/LaAlOs interface remain undetermined. We have mentioned several of these ques- 
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tions in the course of this progress report. In this section, we discuss them in more detail 
and suggest several experiments and computations that may help to further illuminate the 
underlying physics of the range of intriguing behavior observed in this remarkable system. 



A. Whence insulating ra-type interfaces? 

It is found in exp eriment that for samples grown at very hi gh oxygen partial pressure 
(> 10~ 2 mbar) 40j or annealed at 300 mbar during cooling [41] . the n-type interfaces 
can become completely insulating. In fact, the conducting properties of n-type interfaces 
become poorer with increasing po 2 during the growth process (from 10~ 6 to 10~ 3 mbar) 



4=11 ] . With presumably better stoichiometry control, one would expect that the transport 
properties would converge to those of an ideal interface, where the polar catastrophe plays 
the dominant role. What is noticeable is that the high-pressure growth experimental results 
are at variance with these expectations. 

At first glance, this phenomenon would imply that the observed interface conductivity 
comes only from the oxygen vacancies. With the higher po 2 or post-annealing to repair the 
oxygen off-stoichiometry, the concentration of oxygen vacancies should decrease, resulting 
in higher sheet resistance and lower sheet carrier density. However, the sheet carrier density 
does show an abrupt jump at 4 unit cells of LaAlOs [13J, which is difficult to explain using 
oxygen vacancy doping alone. 

We mention here a number of possible scenarios and related questions. The first is that, 
for ideal n-type interfaces and for polar LaAlOs, the polar catastrophe does not take place 
and that the theoretical models are unable to describe the correct physical properties. While 
logically possible, this scenario seems unlikely because the polar catastrophe for idealized 
interfaces is based on simple physical properties such as ionic charges of the cations. Next, 
one can consider the possibility that the theoretical results are correct but inapplicable; that 
the LaAlOs becomes non-polar when the samples are grown under high po 2 or post-annealed. 
Perhaps, aside from removing oxygen vacancies, the oxidation has other side effects. The 
oxidation introduces oxygen atoms at the surfaces of the samples which then diffuse into 
the body of the materials. Does a larger flux of oxygen through the lattice affect the cation 
sub lattice by enhancing disorder and/or intermixing? Does one achieve sharp and ideal 
interfaces or is the disorder sufficient to make the LaAlOs nonpolar? 
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Recent work by Segal et al. on molecular beam epitaxy grown SrTiOs/LaAlOs interfaces 
yield a number of interesting findings concerning some of these questions. Similar to Thiel 
et al. 13], a critical separation of 4 unit cells of LaAlOs is found, as well as similar sheet 
carrier densities in transport. However, x-ray spectra of the La core states do not reveal 
the expected LaAlOs polar field. Quantitatively, the polar field is ~10 times smaller than 
the typical theoretical value of ~0.2 eV/A for samples with LaAlOs thicknesses below and 
above the critical separation. This provides some evidence that the LaAlOs may not be 
polar in some samples. However, it also creates puzzles: if the LaAlOs is not strongly polar, 
then what mechanism leads to a critical thickness close to the expected value? 

One direct experimental probe for these questions is the electron energy loss spectroscopy 
(EELS) performed in the scanning transmission electron microscope (STEM). Detailed stud- 
ies of the Ti, La, and O core level have provided valuable information on the distributions of 
the cations and oxygen vacancies at n-type and p-type interfaces [ijj]. However, we are not 
aware of data published to date on the La core levels for relatively thin LaAlOs films close 
to the critical thickness. The existence of the polar field in LaAlOs will lead to an energy 
shift between neighboring La core levels along the (001) LaAlOs fil m direction. An expected 
shift of about 0.7 eV per LaAlOs unit cell is within present experimental resolution. That 
will allow for direct visualization of the strength of the LaAlOs polar field. 

On the theoretical front, because ideal interfaces and films may not provide a complete ex- 
planation, first principles studies can investigate aspects of disorder and intermixing. Specif- 
ically, which types of cation intermixtures are thermodynamically favored? What effect do 
they have on the polar field (if any)? Do they modify the critical separation or create new 
mechanisms for electron doping independent of the polar catastrophe? In addition, any 
experimental sample will have some degrees of off-stoichiometry. Theoretical studies of the 
possible importance of this disorder effect would be valuable. 

B. What is the value of the critical separation? 

As discussed above, first principles simulations predict an n-type interface critical separa- 
tion of 5 — 6 unit cells of LaAlOs. Recently, based on evidence from optical second harmonic 
generation, Savoia et al. 3] propose that, even at 3 unit cells of LaAlOs, electrons are 
injected into the interface but become localized. The conductance does not appear until the 
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electron gas becomes more uniform at a larger LaA103 thickness. Sing et al. 45J use hard 
x-ray photoelectron spectroscopy and find finite sheet carrier density at samples with 2 u.c. 
of LaA103. They remark that part of the carrier concentration may be of intrinsic origin, 
but it does not lead to conduction. 

What are some possible mechanisms that lead to a smaller critical separation than the 
polar catastrophe would predict? First, oxygen vacancies in the SrTiC>3 substrate are an 
extrinsic doping mechanism that is unlikely here because it can not account simply for the 
abrupt insulating-to-metallic transition. Cation substitution 



which results in the forma- 



tion of metallic Lai-^Sr^TiOa is another possible mechanism. However, some details are not 
clear yet for this mechanism. What is the most stable distribution of cation substitution? 
How does cation substitution depend on the thickness of LaAlOs? In the very thin LaAlOs 
films (< 4 unit cells), is the cation mixing suppressed or enhanced? Cen et al. 21| have 
argued that the creation and annihilation of oxygen vacancies on the surface can induce an 
insulating-to-metallic transition, which raises some new important questions. How large is 
the formation energy of oxygen vacancies on the LaAlOs surface? Does the formation energy 
depend on the thickness of LaAlOs? If thermodynamically favored, does kinetics allow the 
formation of oxygen vacancies on the surface on a reasonable time scale? There are likely 
other possible mechanisms that we have not mentioned here. 

Overall, there appears to an incomplete agreement between experiment and theory con- 
cerning the value of the critical separation. This may suggest that the polar catastrophe 
at the ideal interface, while a reasonable starting point, does not capture all the physical 
effects at play, and that there can be other mechanisms that can generate metal-insulator 
transitions {e.g., oxygen vacancies on the LaAlOs surface). 



C. Is there conductivity on the surface? 

Based on the polar catastrophe scenario, electrons transfer from the surface (or capping 
interface) to the n-type interface when the critical separation is exceeded. In principle, that 
should lead to two regions containing carriers and thus two conducting channels: electrons 
at the n-type interface and holes at the surface (or in the capping layers). However, to date, 
only electron-like carriers have been reported for the n-type interface. This is due in part 
to the fact that it is difficult for simple Hall measurements to detect multi-type carriers. 
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However, Thiel et al. [13[ show in experiment that the conducting layer is not located on 
the LaAlOs surface. 

Although further experiments are needed to confirm the insulating behavior of the sur- 
faces, Thiel's observation gives rise to the following question: if there are holes on the surface, 
why do they not contribute to the transport? One possible answer is self-trapping of the 
holes. The self-trapping can be similar in spirit to that of the self-trapped hole polaron 



proposed by Pentcheva and Pickett 



24| at the p-type interface. Many details remain to be 



clarified. Does the presence of the surface enhance or destabilize self-trapping tendencies? 
What are the effects of surface relaxation? Are relatively large correlation effects (large 
Hubbard U p ) still required for the polaron to form on the surface, and what is the correct 
Hubbard U p 7 

Cen et al. [2l| propose that the presence of oxygen vacancies on the surface can dope the 



interface and make it conducting. In this picture, the electrons associated with the oxygen 
vacancies on the surface are pushed to the interface due to the LaA103 polar field. There 
are no holes on the surface, and the only conducting region is the interface which contains 
the electrons. This prediction is different from that of polar catastrophe mechanism, as the 
defects allow for a new degree of freedom. However, experimentally, it could be challenging 
to distinguish the following two cases: whether there are holes on the surface that do not 
conduct, or whether there are no holes at all. Further experiments may clarify the questions 
regarding holes on the LaAlC>3 surface and could simultaneously give indications of what 
mechanisms come into play for thin LaA103 films. 



D. Are there multiple types of carriers? 



The possibility of multiple types of carriers being 
been suggested in both the published experimental 
Recent experiments also invoke the same picture 



present at the n-tvpe interface has 
191 1 and theoretical 30, y4fl literature. 

rf 

73] . Brinkman et al. [19] speculate that 



the magnetic phenomena they observe are due to spin scattering of itinerant electrons off of 



localized electrons that serve as localized magnetic moments. Seo et al. [73J explain their 
optical measurements based on the simultaneous presence of a low density of high-mobility 
carriers that dominate the transport properties at low temperatures, together with a high 
density of low-mobility carriers localized close to the interface that dominate the optical 
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spectra. The theoretical work uses the Anderson localization picture in 2D to argue that 
the bulk of the carriers, which are spatially localized close to the interface and highly 2D in 
nature, do not contribute to transport due to localization. The localization argument pro- 
vides a possible explanation for why transport measurements on samples grown at higher po 2 
~ 10~ 4 mbar show low sheet carrier densities when compared to the theoretical expectation 
of 0.5e/aQ TO . 

The optical approach of Seo et al. [73J is exciting, as it has the potential to separate 
out different contributions from different types of carriers. To date, these experiments have 
been performed on samples grown at low po 2 ~ 10 -6 mbar, with correspondingly high sheet 
densities of up to 3 x 10 17 cm -2 , beyond the "intrinsic" limit of 0.5e/ag TO . The majority 
of carriers may stem from oxygen vacancies, swamping out any signal from the intrinsic 
electron gas at the interface more evident in samples grown at higher po 2 ■ 

Regarding the Anderson localization scenario, one experiment is to try to grow similar 
samples but to find a way to modify the degree of disorder at the interface. That could 
be induced by intentionally introducing a low density of isovalent cations during growth 
of the interface region or by introducing intentional off-stoichiometry. More quantitatively 
speaking, from the estimation in Section IHIB 31 the required localization behavior requires 
mean free paths on the order of a few nanometers or less, which translates into point-like 
defects every nanometers as well. Since this spacing corresponds to only a few percent areal 
density per interface unit cell, the localization behavior should show strong dependence 
on the density of defects introduced on the percent level in the interface region. On the 
theoretical side, more detailed 2D transport modeling of the n-type interface, including 
scattering and thus localization, can help show whether the picture is applicable and under 
what conditions of disorder. In addition, the expected density and scattering efficacy of 
various defects (e.g., oxygen vacancies or cation intermixtures) could be evaluated from first 
principles simulations. 



V. CONCLUSION 



The SrTiOs/LaA103 interface proves to be an excellent example which shows the richness 
of new phases emerging at complex oxide interfaces. The dependence of each new phase on 
the growth conditions and post-annealing procedures is not yet clear and warrants further 
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study, with a promising outlook towards various device applications of different functionali- 
ties. The SrTiOs/LaAlOs interface paves the way for engineering interfaces between transi- 
tion metal oxides and serves as an archetype of polar-nonpolar interfaces. With advances in 



m 



75j . deepening our 



thin film growth techniques, many new heterointerfaces will emerge 
understanding of interface phenomena and extending the exceptional electronic properties 
in oxides. 
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